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The gapless energy band spectra make the strnctnres based on graphene and graphene bilayers 
with the popnlation inversion created by optical or injection pumping to be promising media for 
the interband terahertz (THz) lasing. However, a strong intraband absorption at THz frequencies 
still poses a challenge for efficient THz lasing. In this paper, we show that in the pumped graphene 
bilayer strnctures, the indirect interband radiative transitions accompanied by scattering of carriers 
cansed by disorder can provide a snbstantial negative contribution to the THz conductivity (together 
with the direct interband transitions). In the graphene bilayer structnres on high-K substrates with 
point charged defects, these transitions almost fully compensate the losses due to the intraband 
(Drude) absorption. We also demonstrate that the indirect interband contribution to the THz 
conductivity in a graphene bilayer with the extended defects (snch as the charged impurity clusters, 
surface corrugation, and nanoholes) can surpass by several times the fundamental limit associated 
with the direct interband transitions and the Drude conductivity. These predictions can affect the 
strategy of the graphene-based THz laser implementation. 


The absence of a band gap in the atomically thin car¬ 
bon structures,such as graphene and graphene bilayers, 
enables their applications in different terahertz (THz) 
and infrared devices One of the most challenging 

and promising problems is the creation of the graphene- 
based THz lasers These lasers are expected to 

operate at room temperature, particularly, in the 6-10 
THz range, where the operation of HI-V quantum cas¬ 
cade lasers is hindered by the optical phonons @ . Recent 
pump-probe spectroscopy experiments confirm the pos¬ 
sibility of the coherent radiation amplification in the op¬ 
tically pumped graphene imii, enabled by a relatively 
long-living interband population inversion 0. As op¬ 
posed to optically pumped graphene lasers, graphene- 
based injection lasers are expected to operate in the con¬ 
tinuous mode, with the interband population inversion 
maintained by the electron and hole injection from the 
n- and p-type contacts [l^. A single pumped graphene 
sheet as the gain medium provides the maximum radia¬ 
tion amplification coefficient corresponding to the quan¬ 
tity dTrcrg/c = na = 2.3%, where ctq = e'^/Ah is the 
universal optical conductivity of a single graphene layer, 
e is the electron charge, h is the Planck constant and 
a ~ 1/137 is the fine-structure constant El- The THz 
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gain in the graphene bilayers [IM or the non-Bernal 
stacked multiple-graphene layers can be enhanced ap¬ 
proximately proportional to the number of the layers. 
However, more crucial is the problem of competition be¬ 
tween the interband radiation amplfication and intra¬ 
band (Drude) radiation absorption [13,113 ■ The latter 
scales with frequency w approximately as 1/(1 -I- w^r^), 
where r is the momentum relaxation time. Hence, the 
onset of terahertz gain is typicall y b elieved to occur only 
in clean samples, where wr ^ 1 |2l| . 

In this letter, we show that the intraband radiation ab¬ 
sorption in disordered graphene bilayers does not actually 
pose a problem for the THz lasing. On the contrary, the 
presence of a certain type of defects reinforces the neg¬ 
ative contribution to the THz conductivity and, hence, 
improves the gain properties of the pumped graphene bi¬ 
layers. The origin of such an effect is associated with the 
indirect interband radiative transitions (see Fig. 1), in 
which the electrons from the conduction band emit the 
photons being scattered on the defects, thus contribut¬ 
ing to the radiation gain. The inclusion of such processes 
can be important in the indirect-gap materials and gap¬ 
less semiconductors. However, in single-layer graphene 
with charged impurities, the indirect interband processes 
contribute to the net THz conductivity a only moder¬ 
ately due to a low density of states (DoS) near the band 
edges (except, possibly, for the case of the carrier interac¬ 
tions with artificial fairly large-size scatterers [13) • The 
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FIG. 1: Energy band diagram of a graphene bilayer with pop¬ 
ulation inversion (left) and schematic energy dependence of 
its density of states (right). Direct interband transitions in¬ 
volve only vertical photon emission (wavy arrow), while in¬ 
direct emission and absorption processes are accompanied by 
carrier scattering (straight horizontal arrows). 


situation is remarkably different is the intrinsic graphene 
bilayers, where the electron-hole dispersion is gapless and 
almost parabolic [l^. With constant DoS near the band 
edges, the indirect interband radiative processes signif¬ 
icantly contribute to the net radiation gain under the 
population inversion conditions. The THz gain is pro¬ 
portional to —Re a. We analyze the spectral dependen¬ 
cies of the real part of the net dynamic THz conductivity 
Re a in the pumped graphene bilayers for the scattering 
by the point ch arg ed impurities (24 . and by the im¬ 
purity clusters . Our analysis shows that in the 

first case, the gain due to the indirect interband transi¬ 
tion almost fully compensates the losses due to the Drude 
absorption. In the case of the cluster scattering, the gain 
he indirect interband transitions exceeds the intraband 
losses. Depending on the cluster size and their density, 
the ratio |Re cr\/<jQ can markedly exceed unity resulting 
elevated net THz gain. 

The real part of the net in-plane dynamic THz con¬ 
ductivity, Re cr, comprises the contributions of the direct 
(vertical) interband electron transitions, Re Ud, and the 
contributions. Re and Re of two types of 

the indirect electron transitions (inside both the bands 
and between the two bands), respectively: 


Re u = Re ad + Re + Re a^T- 


( 1 ) 


The interband conductivity of graphene bilayers due 
to the direct transitions is given by [T^, | 2 ^ 


Re (Jd = aQ ^ [/„(-fej/ 2 ) - fc{hw/2)] 

nw -I- 7i 

hLo-\-2^1 , / hbJ — 2£p\ 

= '^Ql-^^tanh(—^), (2) 


Hui + 


4T 


0.4 eV is the hopping integral between carbon atoms in 
adjacent graphene planes (28| |. Here and in the following, 
for the intrisic graphene bilayers we assume fv and 
are the Fermi functions characterized by the quasi-Fermi 
energies = —ep and /Tc = sp, respectively. This is jus¬ 
tified by strong carrier-carrier scatte ring leading to a fast 
thermalization of excited carriers [a, |l^ . In the THz fre¬ 
quency range of interest, /lo; <C 71 , and the conductivivty 
of bilayer is simply twice as large as that for a single layer. 
Under the pumping conditions, = —fJ-v = and, 

according to Eq. (2), the THz conductivity is negative for 
photon energies huj below the double quasi-Fermi energy 
£ of pumped carriers [i3- 

To evaluate the real part of the THz conductivity due 
to the indirect transitions, we calculate the second-order 
transition amplitudes for the photon emission (absorp¬ 
tion) accompanied by the electron scattering and apply 
the Fermi’s golden rule. As a result, the general expres¬ 
sions for the indirect intra- and interband THz conduc¬ 
tivities read as follows: 


Re ^ E [/A(£k) - /A(ek0] 

k,k^A 

X S{f>w + £k- ek')|Rs(k - k')pnk^(vk' - Vk)^, (3) 


Re ^ T.[fv{e^) - fcie^)] 

k,k' 

X S{hiu - Ek - ek')|Rs(k - k')|^Mkk'(vk' + Vk)^. (4) 

Here A = c, u is the index corresponding to the conduc¬ 
tion and valence bands, Vk = 2hvQ'k./^i is the electron 
velocity in the graphene bilayer, uq — 10 ®cm/s is the ve¬ 
locity characterizing the energy spectra of graphene and 
graphene bilayers. Us(q) is the q-th Fourier component 
of the scattering potential, and = (1 + AA'cos^kk') 
is the overlap between the envelope wave functions in 
graphene bilayers. Equation (3) reproduces the well- 
known result for the dynamic conductivity in the high- 
frequency limit oj ^ T~^, where t is the momentum re¬ 
laxation time, which can be derived from the Boltzmann 
equation. For a correct qualitative description of the low- 
frequency conductivity, we replace the frequency w in the 
denominators of Eqs. (2) and (3) with (w^ -I- 

In the most practical situations,the Coulomb scatter¬ 
ing by the random substrate-induced charged impuri¬ 
ties is the main factor determining the conductivity . 
Considering impurities as random uncorrelated point 
scatterers with the average density rii, we write the scat¬ 
tering matrix element as js^ 


|Rs(q)P = n. 


2'Ke^ 1 ^ 

.K{q + qs)\ 


(5) 


where fvis) and /c(e) are the carrier distribution func¬ 
tions in the valence and conduction bands, and 71 ~ 


where k is the background dielectric constant, and qs is 
the Thomas-Fermi screening wave vector (^ . In the case 
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FIG. 2: Contributions of different radiative transitions to 
dynamic THz conductivity of a graphene bilayer on Si02 
(k = 2.5) and its net value (solid line) versus frequency for 
Hi = 10^^ cm~^ and ef = 50 meV (upper panel), and ef — uj 
color map (bottom panel). Numbers above the lines indicate 
the values of Re (t/ctq. The region above dashed line corre¬ 
sponds to the negative conductivity. 


of the interband population inversion, the latter is found 
to be qs = 4ac7i(l — lii2)//iuo, where Oc = e’^/fiKVo is 
the coupling constant. Evaluating the integral, we find 
the scattering-assisted interband conductivity 


Re (T 


inter 

ind 


STra^CTQ 


vlui 


X tanh 


(w2 + r-2) 
huj — 2eF 


AT 



where $(a:) = dy(l— \1 — y\){l — y/2)/+ x)'^ with 
the following asymptotes: <I>(a:) ~ 21n2 — 1/2, cc <C 1, 
$(a:) ~ l/(2a:2), x :$> 1. At intermediate frequencies, 
l/r <C w <C h{voqs)^/^i, the conductivity due to the 
indirect interband transitions scales as uj~^. This is sub¬ 
stantially different from the case of a single graphene 
layer, where it tends to a constant. A pronounced in¬ 
crease in the indirect interband contribution to the dy¬ 
namic conductivity in a bilayer at low frequencies is at¬ 
tributed to the constant density of states in the vicinity 
of the band edges. 


FIG. 3: The same as in Fig. 2 but for graphene bilayer clad 
by Hf 02 layers (k = 20). 


Figure 2 (upper panel) shows the spectral dependen¬ 
cies of the contributions of different radiative transitions 
to the dynamic THz conductivity as well the net value of 
the latter for a moderate-quality exfoliated graphene bi¬ 
layer [ 2 ^ on the Si 02 substrate with the effective dielec¬ 
tric constant k = 2.5 (graphene bilayer sandwiched be¬ 
tween Si 02 and air) impurity density rii = 10^^ cm“2^ at 
temperature T = 300 K, and for the quasi-Fermi energy 
Ef = 50 meV. The bottom panel in Fig. 2 demonstrates 
the two-dimensional maps of the dynamic THz conduc¬ 
tivity versus the frequency and the quasi-Fermi energy 
for the same structure. As seen, the contribution of the 
indirect interband transitions in such a sample is weaker 
than those of both the indirect intraband and direct in¬ 
terband transitions, although the former is essential for 
compensating, to some extend, the Drude absorption. 
This enables an increase in iRecrj and some widening of 
the frequency range where Re cr < 0. As follows from 
Fig. 2 (bottom panel), the dynamic THz conductivity is 
negative at Ef ^ 40 meV and a;/27r > 10 THz. 

The relative contribution of the indirect interband 
transitions(compared to the ’’normal” intraband conduc¬ 
tivity) increases with an increasing dielectric constant. 
Figure 3 shows the same characteristics as in Fig. 2, for 
a graphene bilayer being clad by the Hf 02 layers, i.e., 
immersed in a media with a fairly high dielectric con¬ 
stant (k = 20). As seen from Fig. 3, the transition to the 
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FIG. 4: Net THz conductivity vesus frequency calculated for 
graphene bilayers on Si 02 with impurity clusters of different 
size Ic- cluster density and charge are ric = 10^^ cm“^ and 
Z = 7, dielectric constant k = 2.5, and quasi-Fermi energy 
£f = 50 meV. Dashed line corresponds to point charged de¬ 
fects with the same charge and density. 


graphene bilayer structures with higher dielectric con¬ 
stant leads to a much stronger contribution of the indi¬ 
rect interband transitions and, hence, to a much larger 
value of |Re a\. The comparison of Figs. 2 and 3 shows 
that the ratio |Re c!\la q for the minimum value of of cr at 
K = 20 is five-six time larger than for k = 2.5. Moreover, 
in the former case Re cr becomes negative starting from 
ep = 20 meV and a;/27r > 5 THz. At higher values of ep 
(i.e., at a stronger pumping). Re cr can be negative from 
the frequencies of about a few THz to a dozen THz. 

The reinforcement of the negative dynamic THz con¬ 
ductivity effect with increasing background dielectric 
constant demonstrated above is interpreted as follows. 
An increase in k results in the reduced Thomas-Fermi 
wave vector q^. This, in turn, leads to switching from 
the strongly screened to almost bare Coulomb scattering. 
From the energy conservation laws it follows that the in¬ 
direct interband electron transitions are favored by a low 
momentum transfer 5, while for the indirect transitions 
within one band q should be large, namely, q > uj/vq as 
the electron velocity in GBL does not exceed vq. Hence, 
for high values of k, the scattering potential behaves ap¬ 
proximately as Vs(q) oc 1/(7, which supports the inter¬ 
band transitions with a low momentum transfer. 

As shown above, the Coulomb scattering by the 
screened point defects cannot lead to the dominance of 
the indirect interband transitions over the indirect intra¬ 
band ones. However, such a dominance can be realized in 


the case of the carrier scattering on sufficiently extended 
scatterers. We consider the charged cluster of size Ic and 
charge as a continuous two-dimensional distribution 
of charge density p(r) = Zeexp{—r‘^jl^)lTrll (since the 
exact density distribution is not very important, for sim¬ 
plicity, we assume it to be Gaussian [23). Solving the 
Poisson equation for this charge density and averaging 
over the random positions of the clusters, we readily find 
the scattering matrix elements [instead of Eq. (5)] 


l^s(q)P = ric 


27rZe^ 


_K{q + qs) 


exp 



(7) 


Here ric is the density of the clusters. In the limit of 
strong screening, g <C Qs, our model of scattering coin¬ 
cides with the widely accepted model of the ’’Gaussian 
correlated disorder” [3ll - l34| with the root-mean-square 

scattering potential \fv^ k, 2^ j nqs- 

Figure 4 shows the spectral characteristics of the dy¬ 
namic THz conductivity in the pumped graphene bilayers 
[ep = 50 meV) with clusters of charged impurities with 
different size Ic- If the scattering by the charged clus¬ 
ters is a dominating mechanism, the absolute value of 
the net dynamic THz conductivity |Re crjcan markedly 
exceed the fundamental ’’direct interband limit” of 2aQ. 
At reasonable values of /c ^ 6 nm and quasi-Fermi energy 
ep = 50 meV, one obtains |Re cr| > 4 ctq. 

In conclusion, we have demonstrated that the graphene 
bilayers with a long-range disorder (impurity clusters 
with reasonable values of the cluster density and size), 
can exhibit a strong negative THz conductivity with the 
span two times or more exceeding the fundamental limit 
2 (Tq. This effect is associated with the indirect interband 
transitions with the photon emission being accompanied 
by the disorder scattering. For the indirect interband 
photon emission to dominate over the ’’normal” Drude 
absorption, the long-wave length Fourier components of 
the scattering potential should prevail. Such kind of scat¬ 
tering potentials can be formed also by extended surface 
corrugations, quantum dots on the graphene bilayer sur¬ 
face @ , and nanohole arrays . 
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